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ABSTRACT 

We simulated both the matter and light (galaxy) distributions in a wedge of the uni- 
verse and calculated the gravitational lensing magnification caused by the mass along 
the line of sight of galaxies and galaxy groups identified in sky surveys. A large volume 
redshift cone containing cold dark matter particles mimics the expected cosmological 
matter distribution in a flat universe with low matter density and a cosmological con- 
stant. We generate a mock galaxy catalogue from the matter distribution and identify 
thousands of galaxy groups in the luminous sky projection. We calculate the expected 
magnification around galaxies and galaxy groups and then the induced QSO-lcns an- 
gular correlation due to magnification bias. This correlation is an observable and can 
be used to estimate the average mass of the lens population and also make cosmolog- 
ical inferences. We also use analytic calculations and various analysis to compare the 
observational results with theoretical expectations for the cross-correlation between 
faint QSOs from the 2dF Survey and nearby galaxies and groups from the APM and 
SDSS EDR. The observed QSO-lens anti-correlations are stronger than the predic- 
tions for the cosmological model used. This suggests that there could be unknown 
systematic errors in the observations and data reduction, or that the model used is 
not adequate. If the observed signal is assumed to be solely due to gravitational lensing 
then the lensing is stronger than expected, due to more massive galactic structures or 
more efficient lensing than simulated. 

Key words: gravitational lensing - method: analytical - methods: numerical - galax- 
ies: clusters: general - large-scale structure of Universe 



1 INTRODUCTION 

It is now very well established that mass concentrations 
on the line of sight of distant sources act as gravitational 
lenses distorting their image. One of the effects is the shear 
of the image and that can be directly observed in drastic 
cases (strong lensing), or statistically when the lensing is 
weak. The other effect is the magnification of sources. Since 
lensing due to a matter overdensity enlarges the solid an- 
gle of the source and conserves its surface brightness, it can 
bring to view sources that would be too faint to observe 
in a magnitude limited survey, but at the same time it di- 
lutes their population density. This phenomenon receives the 
name of magnification bias, and the two competing trends 
it engenders can give origin to both positive correlation or 
anti-correlation between populations of objects with very 
distinct redshift separation. 

Several groups have measured background-foreground 
angular correlation between populations of distant QSOs 
and nearby galaxies or galaxy groups; see Bartelmann & 
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Schneider (2001) for a review, and Guimaraes et al. (2001) 
for a compilation of several results). Galaxy groups are par- 
ticulary interesting because they trace higher density re- 
gions than galaxies in general, and therefore should yield 
a higher cross-correlation signal. Croom & Shanks (1999) 
found a lack of faint QSOs around galaxy groups and in- 
terpreted the anti-correlation signal as due to gravitational 
lensing. Myers et al. (2003) found that the lack of QSOs 
around galaxy groups persists in larger samples, comparing 
QSOs taken from the 2dF QSO Redshift Survey (Croom et 
al. 2004) to galaxies taken from the Automated Plate Mea- 
surement (APM) galaxy survey (Maddox et al. 1990) and 
the Sloan Digital Sky Survey (SDSS) Early Data Release 
(EDR; Stoughton et al. 2002). Using a simple analysis that 
estimates an effective mass for the groups, both Croom & 
Shanks (1999) and Myers et al. (2003) find values that imply 
a high density universe (fl m ^ 1 )• 

Gaztanaga (2003) finds a strong positive cross- 
correlation between bright QSOs and galaxies from the 
SDSS EDR, to which he suggests the interpretation of a 
large anti-bias (b ~ 0.1) on small scales. A similar result 
is put forward by Myers et al. (2005) from the strong anti- 
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correlation found between faint 2dF QSOs and APM and 
SDSS EDR galaxies. 

Both interpretations - high fi m or high anti-bias - 
are allowed by the underlying lensing theory, where the 
background-foreground cross-correlation due to magnifica- 
tion bias depends linearly on the mass density of the uni- 
verse and on an integration over the mass power spectrum. 
If we treat the two possible explanations independently, the 
high f2 m interpretation acts by increasing the overall lens- 
ing weighting factor, whereas the low b interpretation acts 
by increasing the lens mass. Of course both routes are not 
independent, since we expect higher clustering for a denser 
universe. Nevertheless, the results from these works are ex- 
treme, and at face value are in discordance with other ob- 
servations. Although the errors from the referred works on 
the estimate of the universe matter density parameter and 
galaxy bias are large, these results motivate a more elabo- 
rate analysis of the data. 

An analytical theory of the background-foreground cor- 
relation due to weak gravitational lensing was developed by 
Bartelmann (1995), Dolag & Bartelmann (1997), and Sanz 
et al. (1997). Guimaraes et al. (2001) incorporated into the 
formalism more elaborated galaxy biasing and used scale- 
dependent bias defined by the foreground population power 
spectrum. Menard et al. (2003) expanded the weak lens- 
ing approximation to second-order. Jain et al. (2003) used 
the halo model to reformulate the bias dependence in terms 
of the halo occupation properties of the galaxy population. 
Takada & Hamana (2003) used the halo model and NFW 
profile to compute the full non-linear contribution to the 
cosmic magnification statistics. 

The cross-correlation between galaxies or groups with 
QSOs offers a direct opportunity to study and quantify grav- 
itational magnification and the masses of the foreground 
population of lenses. We approached the problem by simu- 
lating both the matter and light (galaxy) distributions, and 
directly calculating for each foreground lens traced by galax- 
ies the angular magnification. We use a large-scale simula- 
tion of the large-scale structure of the universe and a mock 
galaxy catalogue generated from the former to compute the 
gravitational magnification due to the mass associated with 
galaxies and galaxy groups. We use the observational results 
of Myers et al. (2003 & 2005) as study cases, that is, we seek 
to emulate the parameters of these works in our simulations 
and calculations. 

In Section |21 we describe how we generate a simulated 
matter and galaxy distribution, how we identify groups of 
galaxies in the mock catalogue, and test the simulation 
against results from real galaxy surveys. In Section |3] we 
construct magnification maps for the simulated matter dis- 
tribution, and in Section 2] we calculate the magnification 
around galaxy and galaxy groups, and the corresponding 
QSO-lens cross-correlation. In Section|S]we present a simple 
analytic approach to the cross-correlation calculation and 
compare it with the results from the mass-light simulation. 

We show how the cross-correlation can be used to esti- 
mate the average mass of the foreground lenses in Section 
HJ and also compare the mass estimated from the lensing 
results with the mass obtained directly from the matter dis- 
tribution simulation. We discuss our results in Section |7| 




Figure 1. Galaxy redshift distribution. Circles are for the mock 
catalogue and the solid line is for expression ffTl . 



2 MASS AND LIGHT DISTRIBUTIONS 

We generated a simulation of the matter distribution in a 
10 by 75 degree 2 segment of the universe centred at one ob- 
server at zero redshift and extending to z = 1. We use the 
output of the Hubble Volume Simulation (Frenk et al. 2000) 
that has 10 9 particles of mass M part = 2.25 • 10 12 /i _1 M Q 
in a periodic 3000 3 /i _3 Mpc 3 box. We choose the "concor- 
dance model" simulation, which has Qm = 0.3, S1a = 0.7, 
r = 0.21, (78 = 0.90, initial fluctuations generated by CMB- 
FAST, and force resolution of 0.1fr -1 Mpc. 

Our simulation of the matter distribution does not in- 
corporate the evolution of the density fields, but this short- 
coming is not very relevant for our purposes since we are in- 
terested in lensing by structures at small redshift (z < 0.3). 
Structures at larger redshifts, where evolution could be im- 
portant, act mostly as noise for the lensing signal generated 
by the structures at small redshifts, since the large physical 
separation guarantees that they are uncorrelated. One could 
incorporate evolution by using the light cone output of the 
Hubble Volume Simulation, but so far no one has created a 
galaxy mock catalogue from it. 

We generated a mock galaxy catalogue using the sim- 
ulated mass density field and adopting a bias prescription 
for the galaxy population. We used the code of Cole et al. 
(1998), bias model 2, which is a 2-parameter model based 
on the final density field. The mock catalogue generated has 
4 ■ 10 5 galaxies magnitude limited to B < 20.4, mean red- 
shift z = 0.15, and redshift distribution displayed in Figure 
This distribution is well described by the expression for 
galaxy redshift distributions given by Baugh and Efstathiou 
(1993), with (3 = 1.5 and z, = z/1.412: 



N(z) 



/3z 2 



exp 



(1) 



Following Myers et al. (2003), we use the Turner and 
Gott (1976) algorithm to select groups of galaxies from 
our mock galaxy catalogue. The algorithm groups galaxies 
on the basis of their angular density on the plane of the 
sky, ignoring the redshift coordinate. About 45 per cent of 
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Figure 2. Top panel: galaxies (dots) and groups with 7 or more 
members (circles indicate the group centre) in a 60x60 arcmin 2 
sky patch. Bottom panel: convergence map in the same region for 
a source population at z = 1. 



groups with 7 or more members identified by the Turner & 
Gott algorithm have at least 7 members that are physically 
grouped (rather than being chance alignments along the line 
of sight), though 95 per cent of "Turner & Gott groups" 
with 7 or more members trace at least one or more triplets 
of physically-grouped galaxies (Myers 2003). Whether the 
groups identified by the Turner & Gott algorithm are physi- 
cally associated or not, they certainly always represent dense 
projections of galaxies along the line-of-sight, which should 
be hot-spots for lensing. 

The top panel of Figure |2] shows a lxl degree 2 projec- 
tion of the simulated galaxy catalogue and the position of 
groups with 7 or more members identified in it. One can ob- 
serve some correlation between the position of groups and 
regions of high convergence (therefore magnification) in the 
lower panel of Figure [5] The next Section will discuss how 
the lensing map is generated. 

Not all of galaxies that are present in this patch of sky 
shown in Figure [2] of the mock catalogue are visible in the 
projection because many occupy the same position. This 
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Figure 3. Sky density of galaxy groups identified in 2D for mock 
catalogue and APM + SDSS. Filled symbols are for simple mem- 
bership (exactly N galaxies), and open symbols connected by lines 
are for groups with at least N galaxies. 



fact is a consequence of the low mass resolution of the N- 
body simulation. Because the mass particles are very mas- 
sive there are too few of them at small redshifts in relation to 
the number of galaxies. Therefore the galaxy mock code at- 
tributes more than one galaxy to some mass particles, which 
explains why in the top panel of Figure|5|the groups of seven 
or more members and the overall sky patch appear to have 
fewer galaxies than expected. 

The low mass resolution, and the group identification 
algorithm based on friend-of-friends, which identifies neigh- 
boring galaxies inside a given radius, induce a higher iden- 
tification of groups than would be expected. Assigning the 
same coordinates to more than one galaxy may yield a false 
group identification. Figure |3] shows the group density on 
the sky as a function of group membership. The sky density 
for groups with 7 or more members is 6.5deg -2 , which is 
larger than that obtained from real data (APM + SDSS), 
3.8deg~ 2 . This discrepancy may be relevant for the objec- 
tives of this paper and reveals that the galaxies in the mock 
catalogue form more groups than what is observed, even 
though they have comparable sky density (540 deg -2 for 
observed galaxies and 530 deg -2 for mock galaxies). One 
way to compensate for this discrepancy is to assume that 
the real group population is the one with same sky density 
as the observed. Mock groups with 9 or more members have 
this property as compared to observed groups with 7 or more 
members, and we will use those in some of or calculations, 
however the difference is small. 

Figure^Jshows the angular auto-correlation for galaxies 
and groups with 7 or more members, comparing the mock 
catalogue with a combination of the APM and SDSS results. 
Mock galaxies have a lower angular auto-correlation than 
real galaxies for separation less than 1 arcmin (below the 
simulation resolution), and a somewhat higher one for larger 
separation. Groups with 7 or more members identified in 
the mock catalogue have comparable clustering to groups 
detected in the APM survey and SDSS. 
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Figure 4. Angular auto-correlation for galaxies and galaxy 
groups with 7 or more members. Errors are field-to-field. 



3 LENSING MAPS 

The full knowledge of the mass distribution between the 
observer at z = and a source distribution at high redshift 
allows us to calculate the gravitational lensing of sources by 
the intervening mass inhomogeneities. 

In the cosmological context, and assuming the Born ap- 
proximation, the convergence can be calculated by 



k(9) 



W(y)S(0,y)dy 



(2) 



where 8 is the density contrast, y is a comoving distance, 
i/oo is the comoving distance to the horizon, and W(y) is a 
lensing weighting function 



W(y) 



G q (y')fK(y'-y)fK(y) 



a (y) 



Mv') 



dy .(3) 



G q is the source distribution, a is the scale factor, and fx 
is the curvature-dependent radial distance. 

Since our simulation of the mass distribution consists 
of a discrete set of mass points, equation J5J can be more 
suitably written for computational purposes as 



k(9) 



i=l 



S(gj, Zj) 



(4) 



The sum is perfomed over all mass particles in the simula- 
tion falling in the sky patch covered. E(0j, Zi) is the surface 
density of one mass particle in one angular cell of the grid 
adopted and 

2 



D s 



(5) 



4itGD d D ds ' 

where D 3 = o(j/»)/if(y*) is the angular-size distance to 
the source (our maps assume a source surface at z = 1), 
Dd ~ a (yd)fi<(yd) is the angular-size distance to the mass 
particle (deflector), and Dds = a(y 3 )fK(y s — yd) is the 
angular-size distance from the source to the mass particle. 
The empty beam (5(9, y) — —1) value for the convergence 
in expression J5J defines Kmin, which can be also obtained 
from the ensemble average of the terms of equation (g| if 
one uses that (k(9)) = 0. 



The method differs from the multi-plane method in that 
it does not divide the space in cubes and does not project 
the mass particles in planes. Therefore it does not suffer 
from discontinuity problems across space as is the case with 
the multi-plane approach, and also avoids the cumulative 
effective smoothing generated by the multi-plane grids. The 
Born approximation is justified for the lensing systems and 
resolution level that we use here (see Jain, Seljak and White 
2000, and Vale and White 2003). 

The resolution of the convergence map is determined by 
its grid size 6 gr id, and the effective angular resolution of the 
projected mass simulation 9 masa , which can be estimated 
using a similar prescrition used by Menard et al. (2003) to 
calculate the effective smoothing scale 



= / W(y) 



(6) 



D ang (y) 

where W(y) is the normalized weighting function (J^J, l maS s 
is a linear resolution of the mass simulation, and D ang (y) is 
a angular-size distance. Assuming l m ass to be given by the 
force resolution of the Hubble Volume Simulation, we obtain 
Qmass = 0.5 arcmin, which is less then the grid size of the 
convergence map Q gr id = 1 arcmin. 

The convergence map produced has minimum value 
min[fi;(#)] = —0.065, and maximum value max[re(#)] = 0.92, 
which indicates that we are probing regions where depar- 
tures from the weak lensing regime, k <C 1, may be non- 
negligible. 

The magnification map is related to the convergence k 
and shear 7 by 

1 



/i(0) = 



[l- K (6)Y-^(9) 



(7) 



The convergence and shear fields have some common sta- 
tistical properties, for example (7 2 (#)) = (k 2 (0)^, but this 
identity is not valid locally. Although ensemble averages of 
convergence and shear fields cannot be rigorously used inside 
expression Q for the calculation of statistical properties of 
the magnification, one can expand expression J7J in pow- 
ers of k and 7 and them use the ensemble averages of these 
fields to calculate approximately the ensemble average of /x 
and powers of it. This is the approach used, for example, by 
Menard et al. (2003) to calculate magnification correlations 
to second-order approximation. Other approach is to use lo- 
cal approximations, and since we are mostly interested in 
the magnification around identified luminous lenses (galax- 
ies and groups), we use this avenue. This will become clearer 
in the next Section. 

One approximation (Fan and Chiueh 2000) is to assume 
the shear to be negligible in expression (J7J, 



(8) 



In fact, for the magnification generated by a NFW or SIS 
profile the shear will only become important in expression 
£3 very close to the centre, or below than the 1 arcmin scales 
that we can probe in our lensing map given the relative poor 
resolution of our mass simulation on small scales. 

For a SIS profile the shear around the centre of the mass 
distribution is equal to the convergence at the same point 



ksis(0) = JSis(0) = 



2na v Dds 1 
c 2 D a 



(9) 



which implies a magnification 
1 



M = 



1 - 2k\ 



(10) 



More generally, for any circularly symmetric profile one 
could use that the shear is given by 



7(0) = «(< 6) - k(9) 



(11) 



where re(< 6) is the average value of the convergence inside 
a radius 6. 

For our simulated groups with 9 or more galaxies the 
modulus of the average shear in ring around the centre of 
the group is approximately 1/3 of the average convergence 
in the same ring. Therefore for scales that we can probe the 
effective magnification J7J falls between approximations JSJ 
and 110II . being closer to the former. 

The weak lensing approximation, k *C 1, is a first order 
expansion of expression J7J, 



1 + 2k. 



(12) 



We calculated the magnification using various approx- 
imations, but limit to report in our plots the results from 
expressions JHJ and 11211 , which already allows us to see some 
departure from the weak lensing linear regime. Expression 
1121 1 directly yields a null mean field for the magnification 
map since {«;(#)) = 0, but the same does not occur with ex- 
pression (|HJ . It is necessary to renormalize the magnification 
map so that {fi(8)) = 1. One can impose this condition by 
the means of two transformations 



H(0) -» n(0)/p 



p(0)-A*(e)-(£-i). 



(13) 



(14) 



where fi is the mean value of the non-normalized magnifica- 
tion map. For our simulation, and using approximation JSJ, 
we obtained fi = 1.0051, and the difference between the two 
renormalizations is negligible for our results here. 

The importance of considering departures from the 
weak lensing approximation has been highlighted by Barber 
and Taylor (2003), Takada and Hamana (2003), and Menard 
et al. (2003). 



4 MAGNIFICATION AROUND GROUP 
CENTRES AND QSO-GROUP 
CROSS-CORRELATION 

We used our magnification maps and the galaxy groups iden- 
tified in the galaxy mock to calculate the average magnifi- 
cation in annular regions around group centres - see Figure 
|K| One can observe that the weak lensing approximation 
underestimate the average magnification, and that there is 
a strong dependence of the magnification on group mem- 
bership. Error bars are standard deviations from the mean 
value considering total sampling, i.e. all cells from the grid 
are probed. A more realistic error estimate for an obser- 
vational setting would consider a more sparse sampling. A 
rough estimate is to assume that the effective sampling 
used on the simulations is given by the density of cells 
(1 arcmin -2 for our grid), and in a real observation it is 
given by the background source density. This would yield 
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Figure 5. Average magnification around mock galaxy groups. 
The number ranges in the legend are the number of galaxies for 
the sets of groups. Thick lines with filled symbols include depar- 
tures from the weak lensing approximation, and thin lines with 
open symbols are the weak lensing approximation for the mag- 
nification calculation. Errors are the standard deviation of the 



for the study case adopted here, which has a QSO sky den- 
sity PQSO = 43deg~ 2 , a factor yf Pcell/PQSO = 9 larger for 
the magnification error bars. 

The magnification alters the light flux S from sources, 
which implies that if fi > 1 a source that could be too faint to 
be observed in a magnitude-limited survey may be brought 
to view. The number density of sources in the magnitude 
limited survey is N(> S/fj,). But at the same time the mag- 
nification also alters the area behind a lens, A' — A/ fx, so 
in a fi > 1 region the area behind the lens is expanded and 
the source counting density is diluted. These two competing 
effects are compared in the net enhancement factor 



A'N(> S/fi) 
AN{> S) 



(15) 



If the source cumulative number counts by flux is of the 
form N(> S) oc S~ a then the enhancement factor reduces 
to q = /i s_1 . The coefficient s of the number-flux relation 
relates to the coefficient of the number-magnitude relation 
by s — 2.5/9. We use (5 = 0.29 where required (Myers et al 
2003). 

The enhancement can be estimated as the ratio of the 
the observed number of QSO-group pairs, DD(0), to the 
expected number of random pairs, DR(9). If we recall the 
cross-correlation estimator 



DR(9) 



(16) 



then 



u 9g (d) = Key 1 - 1 ■ (i7) 

Figure |5] shows the QSO-group correlation function for 
groups with different galaxy membership. A departure from 
weak gravitational lensing is illustrated for groups with more 
than 20 members. The curves show a strong dependence in 
relation to group membership, indicating its relation to the 
underlying mass overdensity. 
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Figure 6. Cross-correlation between QSO and mock galaxy 
groups (we plot —u qg to allow the use of logarithmic scale). The 
number ranges in the legend are the number of galaxies for the sets 
of groups. Thick lines include departures from the weak lensing 
approximation, and thin lines are the weak lensing approximation 
for the magnification calculation. Curves for groups of member- 
ship 1 to 4 are shown individually (from bottom to top for low 
to high membership). Errors are the standard deviation of the 
mean. 



Figure |7| compares the cross-correlation results from 
simulations to data from Myers et al. (2003) for groups, 
and Myers et al. (2005) for galaxies. Simulation results for 
angles smaller than 1 arcmin cannot be obtained due to 
limited simulation resolution; however for angles from 1 to 
100 arcmin the comparison with data gives a large dis- 
agreement between the amplitudes of observed and simu- 
lated cross-correlations. We find the parameters of a power 
law that best describe the observed and simulated QSO- 
galaxy and QSO-group cross-correlations, and quantify the 
disagreement in amplitude between the observed and simu- 
lated results, which is larger for groups (factor of ~20) than 
for galaxies (factor of ~7). 



5 SIMPLE ANALYTIC APPROACH 

In this section we present a simple analytic approach to 
check and help the interpretation of the results from sim- 
ulation. 

Expression llh^ for the cross-correlation is equivalent to 
(see Appendix [XJ 



n q ((f>) 



- 1 



n g (<t> + 9) 



- 1 



(18) 



where n q and n g are the QSO and galaxy (or galaxy group) 
densities (a bar over a quantity indicates its mean value), 
and (...) represents the average over </> and the direction of 
9 (but not its modulus). 

From expression I18H , and using the formalism of Dolag 
and Bartelmann (1997), one can derive assuming weak lens- 
ing that 
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Figure 7. QSO-galaxy and QSO-group cross-correlation. Obser- 
vational data is from Myers et al. 2003 and 2005 with field-to-field 
errors. Simulation uses groups with 9 or more galaxies and esti- 
mated error with same source density as observed data. Some 
observational points fall below the shown logarithmic scale 



x/ dkkP gm (k,y) J [f K (y)k6] , 



(19) 



where y is the comoving distance, which here parameter- 
izes time (tjoa represents a redshift z = oo), and k is the 
wavenumber of the density contrast in a plane wave ex- 
pansion; Jo is the zeroth-order Bessel function of first kind; 
and }k (y) is the curvature-dependent radial distance (= y 
for a flat universe). P gm (k,y) can be seen as the galaxy- 
mass cross-power spectrum (Jain et al. 2003), and under 
some assumptions (Guimaraes et al. 2001) may be expressed 
as P grn (k,y) = y / P g (k)P m (k,y), where P g (k) is the power 
spectrum for galaxies or galaxy groups and P m (k,y) is the 
non-linear time evolved mass power spectrum. 

Menard et al. (2003) calculated the second order con- 
tribution to ui qg (8) in addition to the dominant first order 
term given by expression 11911 and found that it can increase 
the amplitude of ui qg (9) by 15% to 20% on scales below one 
degree. 

The mass power spectrum P m (k,y) for the ACDM 
model can be obtained analytically from the linear theory 
and the use of the non-linear prescription given by Peacock 
& Dodds (1996). Figure |H] shows the mass power spectrum 
at 2 = obtained by the method described. 

To obtain the three-dimensional power spectrum in real 
space for galaxies and groups P g (k) we use the corresponding 
two-point angular auto-correlation function iu gg (9) and the 
relation (Peacock 1991) 



dy y 4 4> 2 {y) 



dk^-P a {k)J {ky6) , (20) 
Zty 



where 4>{y) is the selection function, normalized such that 
J y 2 <j>dy — 1. If we assume (f> oc y 1 ^ 2 exp[— (y/y*) 2 ], uJ gg (9) = 
B8 13 , and that the power spectrum is a power law, then 
expression (1201 can be inverted, allowing us to find that 



P g (k) = 



■^yl~ 



B 



(/? + 2)r(-f)r»(f) 



-(3+2) 



(21) 
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Figure 8. Power spectra for galaxies, groups, and mass distribu- 
tion. The mass power spectrum is for a concordance model. The 
power law spectra for the APM+SDSS and mock catalogue sets 
are obtained from the respective auto-correlation function using 
the method described in the text. 



We find that y, = 400/i -1 Mpc allows a good approxi- 
mation for the selection function used in the generation of 
the galaxy mock catalogue. 

We plot on Figure |H1 the power law spectra calculated 
from the auto-correlation functions shown in Figure 01 using 
equation 1211 . and also the original spectra for APM galaxies 
(Gaztanaga & Baugh 1998) and Abell-ACO clusters (Miller 
& Batuski 2001) for comparison. 

Figure [§] shows the QSO-galaxy and QSO-group cross- 
correlation obtained from using Equation 1191 and the ap- 
proximated power spectra for galaxies and groups obtained 
from the auto-correlation functions, using the method de- 
scribed in this section. There is a general agreement between 
our mass-light simulation method results and the simple an- 
alytic cross-correlation results, which is a useful confirma- 
tion of the consistency of both approaches. 

The analytic cross-correlation expression 1191 can be 
further simplified if we approximate the source and lens dis- 
tributions by Dirac Delta Functions peaked at the average 
comoving distances (y q for the QSOs and y g for the galax- 
ies), and the power spectrum term by P gm (k) = Ak a , then 



w 9ff (0) 



3-2 a (s-l) (H \ 2 F(f + 1) 



Jk (y g ) - fK(y g ) 



a (y s )fK(y q )fK (y g 



r -4 



-n m A0- 



(22) 



where 6 is in radians. If one assumes linear bias, then from 
expressions 1221 and 1211 it can be deduced that the cross- 
correlation has the same angular dependence as the galaxy 
or group auto-correlation, that is to say, if uj gg (9) oc 9 13 then 



9 1 . This simple relation is roughly confirmed by 



the results (both from simulation and observation) shown 
in Figure 0] for the auto-correlation for galaxies and groups 
and in Figure |7| for the respective cross-correlations. 
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Figure 9. QSO-galaxy and QSO-group cross-correlation from 
simulation and analytic calculation. For the analytic results the 
only difference in the calculation is the lens (galaxy or group) 
power spectrum used. 



6 MASS ESTIMATION 

One simple way to estimate the mass of a lens is to assume a 
mass profile for it, calculate the expected magnification, and 
use expression 1171 to determine the cross-correlation that a 
population of these lenses would generate (halo lens fitting). 
One can then find the parameters for the mass profile chosen 
that best fit the observed QSO-lens cross-correlation, and 
with it the mass comprised by the lens. 

This lens mass estimation method was used by My- 
ers et al. (2003) in groups of galaxies identified in a two- 
dimensional sky projection, however these groups are not 
necessarily bound entities. Also, the method essentially as- 
sumes that the lenses are identical and isolated objects, ig- 
noring large-scale structure considerations such as cluster- 
ing and filaments. Nevertheless, the simplicity of the model 
is appealing, and its result can be attributed to some "ef- 
fective" halo. Our simulation allows us to test the validity 
of the method by comparing our mass estimates for lenses 
obtained from the QSO-lens cross-correlation and directly 
from the mock matter distribution. 

We examine two popular choices for halo profiles, 
the singular isothermal sphere (SIS) and the NFW profile 
(Navarro, Frenk and White 1997). 

For a SIS halo profile psis(r) = rr v /2-ivGr 2 , and the 
mass inside a radius r is Msis = 2a 2 r/G, where a v is the 
velocity dispersion. 

For the NFW halo profile 

pNF\v{r) = , (23) 

(r/r s )(l + r/r s ) z 

where p c is the critical density. We reduce the NFW profile 
to a one parameter description, the mass inside a 1.5/i _1 Mpc 
radius sphere (Mi. 5), using the relation (Maoz et al. 1997) 



0.3 



M1.5 
10 15 M C 



/i _1 Mpc 



(24) 



where we use A = 1/3. 

The convergence k(9) generated by a halo can be cal- 
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Figure 10. Excess surface density Tj exc around galaxies and 
groups in relation to the mean surface density of the universe 
to redshift 0.4 (4.3 • 10 12 ft _1 M e arcmin 1 ). The legend numbers 
indicate the galaxy group membership. Errors are the standard 
deviation of the mean. 

culated projecting its mass density profile into a plane and 
using equation Q. The shear can be obtained in this case of 
circular symmetry by equation I1H . See Wright & Brainerd 
(2000) for explicit analytical expressions for lensing by ha- 
los. The magnification can then be calculated using equation 

For the galaxies and galaxy groups in our mock cat- 
alogue we can estimate the lensing mass associated with 
them directly from the matter simulation. Figure ITTT1 shows 
the excess surface density E elc in a ring around galaxies and 
groups in relation to the mean surface density of the universe 
to redshift 0.4. The excess density is largely insensitive to 
redshift cuts above this value, since structures above it are 
not correlated to the luminous structures that we are inter- 
ested. 

The excess mass contained inside a redshift cone centred 
at a lens (obtained from an integration of the excess surface 
density) can be used to give an estimate of the average lens 
mass. However this method in fact includes mass outside 
the nominal radius, since the cone angular aperture only 
excludes the mass outside the radius in the perpendicular 
direction of its axis. The projected mass outside the radius 
in the line of sight direction is also included. Therefore one 
would expect the mass obtained through this method to be 
an over-estimate of the real mass inside a chosen radius. We 
show these masses around lens centres for cones with angular 
aperture corresponding to 1.5/i _1 Mpc radii for galaxies at 
an average redshift of 0.15 in Tabic Q 

Table compares the masses associated with galaxies 
and groups using halo lens fitting to the simulation and ob- 
servational QSO-group cross-correlation results and directly 
from the matter distribution simulation. 

For galaxies the redshift is well defined (to projected 
galaxy groups it is not) and therefore we can estimate the 
mass inside a sphere centred on the galaxy coordinates. We 
calculate this mass for a 1.5/i _1 Mpc radius from the cen- 
tre of the lenses. Such a large radius for a galaxy is chosen 
to make possible the comparison with similar measures for 



groups and therefore should be seen as a measure not of the 
individual galaxies but of the galactic environment. 

Looking across the row in Table Q for the simulation 
results having galaxies as lens centres, we note that the av- 
erage mass inside a sphere of radius 1.5/i -1 Mpc is roughly 
six times smaller than the mass computed from the excess 
surface density. The masses estimated from the simple halo 
profile lensing models fall between these two. Assuming that 
the same relation between the described masses applies for 
groups, then the estimated masses from the halo lens fittings 
are in reasonable agreement (tending to be an overestimate) 
with the average real mass associated with groups. 

The vertical comparison of the values in the three first 
columns in Table^also reflects the significant amplitude dis- 
agreement discussed in Section [I] between observation and 
simulation, but to a less dramatic degree than the power law 
analysis. For the galaxy results the mass inferred from ob- 
servations using the halo profile lensing method is ~5 times 
larger than the corresponding results from the simulation, 
and for groups this same factor is ~10. 



7 DISCUSSION 

We have pursued a computational approach to the prob- 
lem of QSO-group cross-correlation, which uses the same 
galaxy group identification procedure of Myers et al. (2003) 
in galaxy survey data, who found a large anti-correlation sig- 
nal. The method consists in simulating the mass and light 
distribution, and directly computing the magnification due 
to gravitational lensing caused by the mass concentrations 
traced by selected galaxies. The background-foreground cor- 
relation is computed directly from the average angular mag- 
nification around the group centres, while previous works 
that also use simulations calculate the cross-correlation from 
lensing field statistics. Our approach takes into account the 
large scale structure of the universe, deviations from the 
weak lensing approximation, the actual galaxy distribution, 
and selection criteria. 

The large scale structure of the universe implies that 
the groups are clustered, and therefore this feature is prop- 
erly taken into account in our calculations, in contrast to 
a simpler approach that models lenses as isolated objects, 
as used by Myers et al. (2003). Nevertheless, this simple 
approach was shown to be roughly adequate for lens mass 
estimation. The mass obtained from isolated halo fitting to 
cross-correlation results obtained using mock galaxies drawn 
from a suitably biased Hubble Volume simulation are com- 
parable for the two halo profiles examined, and the direct 
mass estimate from the mock matter distribution. 

The consideration of the LSS also allows the inclusion of 
structures that are not connected to the groups, but act as 
noise generators for lensing measurements. In fact, most of 
the lensing for distant sources come from structures around 
half of the way to the observer, and these structures are not 
connected to the galaxy groups since they are much closer 
to the observer. This background lensing signal from uncor- 
rected structures is the main noise source for the lensing 
signal from structures associated with observable galaxies. 

Our results for galaxy groups predict a strong depen- 
dence of the cross-correlation on the group galaxy member- 
ship, which is in accord with the calculated mass estimates 
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Table 1. Estimated average masses within 1.5h~ x Mpc of the galaxy or group centre in units of 10 14 /i — 'Mq. The columns are: cr v is 
the velocity dispersion for a SIS model and Mi. 5 sis is the corresponding estimated mass, Mi 5 nfw is the estimated mass in a NFW 
profile model, Mi 5 s exc is the estimated mass from an integrated excess surface mass density around galaxies or groups at the average 
redshift z=0.15 (errors are standard deviation from the mean), and Mi .s jSp h is the average mass inside spheres centred at galaxy positions 
(errors are standard deviation). Mi.5^ exc and Mi .5, sp h are clearly only available for simulations. Groups from simulations have 9 or 
more galaxies and groups from observation have 7 or more members. The observational data sets come from Myers et al. (2003 & 2005). 



data set 


a v (km s - 1 ) 


Miosis 


Mi^NFW 


Mi.s^Ecxc 


Ml.5 iSp h 


galaxy, simulation 
galaxy, observation 


227±f s 
470±g 


0.4±0.1 
1.5±0.4 


n 9 o+0.20 
u - zo -0.19 

1.3±0.5 


0.64±0.01 


0.11±0.10 


group, simulation 
group, observation 


409+« 


1.1±0.3 


l 1+0.6 

12.±9 


1.5±0.1 





for these groups. The cross-correlation curves suggest effec- 
tive mass profiles for the groups that are distinct in ampli- 
tude, but comparable in slope, just favouring slightly cuspier 
profiles for more massive structures. 

The results for galaxy groups of membership 1, i.e. iso- 
lated (field) galaxies, when compared to the results of the 
average galaxy suggest an interesting application. Both the 
QSO-galaxy cross-correlations and the direct mass estimates 
for the two populations of lenses indicates that the average 
galaxy is twice as massive as field galaxies. If the galaxy pop- 
ulation is large enough, then it could in practice be broken 
down into smaller sub-sets according to a chosen character- 
istic and the cross-correlation observation could be used to 
compare the relative masses of these sub-populations. 

Our simulation lacks mass resolution. A higher resolu- 
tion simulation of the matter distribution at small redshifts 
would allow a more precise determination of deviations from 
linearity of the weak lensing approximation. It would also 
allow the investigation of smaller angular scales, including 
halo substructure and the magnification closer to the core 
of galaxies and clusters. It is possible that substructure may 
make non-linear magnification important even at regions far 
from the identified galaxy or galaxy group centre. Therefore 
a higher mass resolution simulation could yield a stronger 
expected cross-correlation due to non-linear lensing, not only 
for small 6 but also for larger angles. This possible effect is 
worth further investigation. 

The persisting discrepancy between the amplitude of 
the expected QSO-galaxy and QSO-group cross-correlations 
in relation to the much lower values expected from the cos- 
mological and lensing models adopted is a stimulus for fur- 
ther investigations. 

It may be possible that there are observational issues 
that cause a systematic error in the cross-correlation de- 
termination. The hypothesis of dust is by now very weak- 
ened, since the reddening of lensed QSO was found not to 
be significant by Myers et al. (2003), and the strong posi- 
tive cross-correlation measured by some groups when using 
bright QSOs also argues against it. The dust would need 
to obscure only faint QSOs and cause no reddening to be a 
viable explanation. 

On the other hand, if all the observational issues are 
very well understood, and the cross-correlation signal due 
to lensing is stronger than predictions such as we provide 
in this paper, then various possibilities are indicated. One 
possibility is a higher mass density for the universe, which 



would increase the overall weighting of the lenses (the lens- 
ing factor). This possibility is highly constrained from other, 
probably more precise and accurate determinations of Q m - 
Another possibility is that the lenses themselves (galaxies 
and clusters) are more massive than the ACDM model pre- 
dicts. A third option is that the lensing efficiency is higher 
than expected from our simulation or analytical calculations 
that assume weak lensing due to a more prominent role of 
non-linear lensing. At high density regions (near the core 
or due to substructure) the magnification is very non-linear 
and could in principle give a higher contribution to the signal 
than what is currently being simulated. This last possibility 
is favored by our observation that the disagreement is larger 
for galaxy groups than for individual galaxies. Since groups 
trace regions of higher density than isolated galaxies, it is 
expected that non-linear effects would be more pronounced 
for groups. 
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APPENDIX A: CROSS-CORRELATION DEFINITION AND ESTIMATOR 

Here we show the equivalence between expressions 116H and 11811 for the angular cross-correlation between elements belonging 
to two populations. 

We can describe the populations of objects by two sets of position vectors a, (i = 1, N a ) and (3- (j — 1, Np), where N a> p 
is the number of elements in either set. So the population density n a can be written as 

N a 

n a (<j>) =^2,8 D {4> - a.i) , (Al) 
where Sd is a Dirac delta function, and the population mean suface density n a is 

n a = (Tt n ((ft)) = j.^ — = — , (A2) 

where A is the survey area, and a similar expression holds for np. 

Therefore the cross-correlation between populations a and /3 is, according to expression 1181 . 

i E«M*-l/3 3 --«< 



n a np 2n6A 
DD{9) 

^ 1 



1-1 + 1 



DR{6) 
which is expression 11611 . 

The last step takes two facts into consideration; that DD(9) is the number of pairs of elements from the two populations 
such that their separation |/3 ■ — on\ is 8, and that for a random distribution of both populations DR(9) = 2n6An a np. 



